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We present Monte Carlo simulations of the formation 
and structure of polysilane, hybrid polysilane/polysilyne and 
polysilyne networks. The simulation technique allows for the 
investigation of large networks, containing up to 1000 silicon 
atoms. Our results show that ring formation is an impor- 
tant factor for all three types of materials. For polysilyne 
networks, a random structure is found incorporating cyclic 
substructures, linear chains and branching points. 



I. INTRODUCTION 

The electronic and optical properties of silicon-based 
materials are strongly related to the structure of the 
silicon skeleton and the size and surface properties of 
the materials O. Crystalline silicon, which consists 
of a three-dimensional silicon framework, is an indirect 
semi-conductor with a band gap of around 1.10 eV. In 
other forms of three-dimensional silicon, such as nano- 
crystalline and porous silicon the band gap is somewhat 
larger and quasi- indirect, so that rather efficient photo- 
and electroluminescence have been observed. In con- 
trast, one-dimensional polysilanes (-SiR2-) n , which are 
linear polymers consisting of a backbone of silicon atoms 
to which two organic side groups R are bonded ^,|| ex- 
hibit well-defined, intense absorption and emission bands 
in the UV spectral region. An interesting aspect of the 
electronic spectra of polysilanes is that both in the solid 
state and in solution they are often strongly tempera- 
ture dependent. This thermochromism finds its origin in 
conformational changes of the silicon backbone, which is 
said to be er-conjugated. Polysilanes are usually synthe- 
sized by polymerization of SiR2Cl2 monomers, employ- 
ing alkali metals as coupling agent. The heterogeneous 
polymerization process is rather complicated, but is pro- 
posed to involve both silyl radicals and silyl anions as 
reactive intermediates Q]. It is noteworthy that in the 
alkali metal-mediated polymerization there is a tendency 
to form cyclic oligosilane oligomers. 

Organic side chain appended silicon-based materials 
with a higher dimensionality than polysilanes can be ob- 



tained by polymerization of trichlorosilanes RSiCl3 ||,[| . 
This yields so-called polysilynes, in which each silicon 
atom is bonded to three other silicon atoms and to one 
organic side group R. Polysilynes exhibit a broad indirect 
semiconductor-like absorption which tails into the visible 
region ^-^|. The fluorescence, emanating from trapped 
excitons, is situated in the visible region and has a broad 
appearance. Hybrid polysilane/polysilyne networks, ob- 
tained by polymerization of a mixture of dichloro- and 
trichlorosilanes, have also been reported |s|-|l^]. These 
kind of polymers can be viewed as linear polysilanes in 
which SiR branching points are incorporated. In hybrid 
polysilane/polysilyne networks the photophysical proper- 
ties of polysilanes and polysilynes are more or less com- 
bined. 

Since their discovery, the structure of polysilynes and 
hybrid polysilane/polysilyne networks has been the sub- 
ject of debate. Polysilynes were initially assumed to con- 
sist of a rigid network of interconnected ring-like struc- 
tures H , but it was also argued that they form essentially 
two-dimensional sheetlike networks . According to an- 
other point of view it was reasoned that the growth of a 
polysilyne preferably occurs at the termini of the polymer 
chain, which leads to a hyperbranched, dendritic mor- 
phology fjlj. Cleij et al. used PM3 calculations on silyl 
radicals and silyl anions to show that chain propagation 
preferentially occurs at the termini of oligo- and poly- 
mers and is more likely to occur than branching [fl4|| . 
It furthermore appeared that formation of silicon ring 
structures is also an aspect to take into account. Hence, 
the PM3 calculations led to the idea that polysilynes pos- 
sess a predominant one-dimensional structure with small 
branches and incorporated rings. 

A number of experimental observations also indicates 
that polysilynes and hybrid polysilanes/polysilynes can 
be viewed as essentially linear structures. For a series 
of n-hexyl substituted hybrid polysilyne/polysilane net- 
works thermochromism (both in the solid state and in 
solution), fluorescence and a degree of exciton delocal- 
ization were observed which resembled the properties of 
linear polysilanes j^] . Even more surprising, for a polysi- 
lyne with oligo (oxyethylene) side chains in aqueous envi- 



1 



ronment thermoresponsive behavior very similar to that 



of linear polysilanes was found 14 1. This thermorespon- 
sive behavior has to originate from folding and unfold- 
ing processes of the silicon framework. These results im- 
ply that the silicon backbone of polysilynes and hybrid 
polysilyne/polysilanes are to a certain extent flexible and 
behave much like one-dimensional systems. Similar con- 
clusions can be drawn from the photophysical properties 
of well-defined oligosilane dendrimers (l^Jll| ■ 

In this contribution, we present computer simula- 
tions of the formation of polysilane, hybrid polysilane- 
polysilyne networks and polysilynes, and study the prop- 
erties of the resulting networks. There are three es- 
sential characteristics that we included in these simula- 
tions: 1) the starting point is a random mixture without 
any polymerization; 2) monomers diffuse, and when they 
meet they can form stable bonds; 3) conglomerations of 
bonded monomers are not stationary and rigid, but they 
show some diffusion and flexibility. More specifically, in 
our simulations we implemented the time evolution by 
means of Monte Carlo dynamics (since the dynamics is 
overdamped because of the solvent), and the elastic prop- 
erties of conglomerations were described by an empir- 
ical interaction potential featuring bond-stretching and 
bond-bending, with parameters chosen to match experi- 
mental properties of crystalline Si. It is anticipated that 
a thorough knowledge of the structure of polysilanes and 
polysilynes gives more insight in the properties of these 
silicon-based materials. 
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where the summation runs over all pairs and triples in 
the system; a and (3 are the bond-stretching and bond- 
bending force constants, respectively; is the vector 
pointing from particle i to particle j and is the 
ground-state distance between particles i and j; if these 
are both silicon atoms, dij equals r$s, otherwise one par- 
ticle is a silicon atom and the other an alkyl fragment, 
and dij equals the equilibrium Si-R distance r^s m that 
case. The values for the parameters used in the simula- 
tion are listed in Table 0. 






II. MODEL OF POLYSILANE AND POLYSILYNE 

The simulations are started without any trace of poly- 
mer present. The only bonds initially present in the sys- 
tem are therefore the internal Si-R bonds shown in Fig. ^. 
However, in the course of the simulation, as monomers 
react with each other, Si-Si bonds also form. In these 
simulations, both the silicon atoms and the alkyl groups 
are treated as hard spheres with radii r$ and t\r, respec- 
tively, see Fig. [|. In the ground state, the Si-R distance is 
set to trs and the R-Si-R bond angle to the tetrahedral 
angle 8o, with cos(0o) = —1/3. During the simulation, 
as the monomers diffuse and react, distances and angles 
are allowed to fluctuate around their ground state values. 

In the model, all non-bonded particles interact accord- 
ing to hard-sphere potentials. Bonded particles interact 
according to the Keating potential |17fl , which describes 
the stiffness of the network with respect to bond-length 
and bond-angle distortions. This potential requires an 
explicit list of all bonds, and is given by 
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FIG. 1. Diagrammatic representation of S1R2 (top) and 
SiR fragments (bottom). Silicon atoms are labeled 'Si', alkyl 
groups 'R'. Their respective radii are rs and tr. The ground 
state Si-R distance is set trs and the ground state R-Si-R 
angle in silane to the tetrahedral angle Qq. 



III. SIMULATION DYNAMICS 

The simulation starts with a configuration of 712 SiR2 
and 713 SiR monomers, placed at random positions and 
with random orientations in a periodic box, under the 
constraint of the hard-sphere interactions. The size of 
the box is chosen such that the particle density equals 
2.0 • 10~ 4 monomers per A 3 , matching the experimental 
conditions where typically 20 mmol of monomer is put 
into a reaction volume of 60 ml ||. The nature of the 
dynamics is two-fold: monomer diffusion and reactions 
between monomers. We simulate this with the approach 
of rare-event dynamics Jjl|. In this approach, the two 
types of events are monomer displacements, which occur 
with a rate r\, or bond formation between two silicon 
atoms, which occur with a rate r^,. This reaction can 
only occur when (1) both Si atoms are not saturated, 
i.e., they have less than four bonds to either alkyl-groups 
or other silicon atoms, (2) the silicon atoms are separated 
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less than a cut-off distance r c , and (3) the silicon atoms 
are not mutually bonded already. We continue to simu- 
late until reactions become rare. This occurs when ap- 
proximately 80% of all possible Si-Si bonds have formed. 
We then investigate the structure of the final network. 

When concentrating on the final structure, the overall 
scaling of the rates r\ and T2 is irrelevant, since it only 
affects the time scale of the simulation. The networks 
presented here were generated using T\ : 7*2 = 5 : 1. To 
the best of our knowledge, no experimental data on the 
reaction rates r\ and ri are reported. However, we have 
established that the networks are insensitive to the ratio 
of rates. For r\ : r2 — 1 : 2 and 1 : 10, we obtained 
similar results. 

We recall that polysilanes are formed by silicon atoms 
bonded to two alkyl groups, and polysilynes by silicon 
atoms bonded to one alkyl group. Therefore, a system 
consisting of n 2 S1R2 and 113 SiR monomers is represented 
by a total number of N = 3ri2 + 2ns locations in three- 
dimensional space. We denote the total number of pairs 
of silicon atoms that are able to react as P, a quantity 
that varies during the simulation. In terms of N and P, 
the total diffusion rate i?i can be written as R\ = r\N, 
the total reaction rate R2 as R2 = T%P, and the total rate 
of events as R\ + R2. Events are selected one-at-a-time. 
To each event, a time-increment of At = l/(i?i + R2) is 
attributed. The likelihood that this event is a diffusion 
event or a reaction event, is equal to R±At and ii^Af, 
respectively. 

To describe the diffusion process, we introduce the fol- 
lowing event: 

1. We select randomly one silicon atom or alkyl group. 

2. For this silicon atom or alkyl group, a displacement 
is proposed, drawn randomly from within a sphere 
with radius r m . 
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If the hard-sphere constraints are violated, the pro- 
posed displacement is rejected. Otherwise, the dis- 
placement is accepted with the Metropolis accep- 
tance probability p9(j 
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where kb is the Boltzmann constant, T is the tem- 
perature, and Eb and Et are the total (Keating) 
energies of the system before and after the random 
displacement. 

Since we displace single silicon atoms and alkyl groups, 
the above procedure also allows for the vibration and 
rotation of monomers. 

To describe reactions in the simulation, we introduce 
the following reaction event: 

1. We select randomly one of the P pairs of silicon 
atoms able to react. 



2. A bond is placed between the two silicon atoms 
constituting the selected pair. 



IV. RESULTS 

Polysilane networks, polysilyne networks and hybrid 
polysilane-polysilyne networks are obtained using the 
simulation method described above. Each system con- 
tains a total of 1000 monomers. To characterize the 
structure of these networks, we proceed as follows. First, 
we identify clusters, defined as a group of connected 
monomers. For each cluster, we calculate its genus g 
given by g = 1 + e — n, where e is the number of Si-Si 
bonds in the cluster and n the number of silicon atoms in 
the cluster. The genus helps to identify the topology of 
the cluster: it measures the number of bonds that can be 
cut before the cluster loses its connectivity. For example, 
if g = 0, the cluster is a chain of connected monomers, 
i.e. a polymer, and if g = 1, the cluster is a ring, possibly 
with a number of side-chains attached. Networks with 
high values of g have a complicated topology. 



A. Polysilane 

We performed one simulation containing 1000 SiR.2 and 
no SiR fragments at a density of 2.0 • 10 -4 monomers per 
A 3 . The final network is shown in Fig. ^; 86.9% of all 
possible Si-Si bonds are formed. 

TABLE I. Values of parameters used in the simulation: 
vrs and rss are the ground state Si-R and Si-Si distances, 
respectively, and rs and tr are the radii of a silicon atom 
and an alkyl group, respectively; values for these four pa- 
rameters were adapted from an MM2-calculated structure of 
dodecamethycyclohexasilane. Reactions between two silicon 
atoms can occur only when their separation is less than r c . 
r m is the maximum particle displacement during a diffusion 
step, a and f3 are the standard Keating parameters jl7j. T 
is the temperature; 7*1 and ?~2 are the rates of diffusion and 
reaction, respectively. 
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Since silane is bifunctional, only two types of structure 
exist in this case: chains (polymers) with g = and 
rings with g = 1. The network consists of 228 separate 
clusters, ranging in size from one up to twelve silicon 
atoms. Of these clusters, 131 are chains, the remainder 
are rings. Polymer length and ring size statistics are 
shown in Fig. ||. In our simulation, 46.8% of the silane 
monomers form linear structures; a substantial fraction of 
53.2% of the monomers is found in rings. For this system, 
the Si-R and the Si-Si bond lengths are 2. 10 ±0.07 A and 
2.32 ± 0.08 A. The R-Si-R and Si-Si-Si bond angles are 
107.0 ±5.1 and 109.7 ± 7.9 degrees, respectively. 




the system contains 104 separate clusters, ranging in size 
from two up to 36 silicon atoms. Of these clusters, 77 
are incorporated in linear structures, the remainder form 
rings. Fig. ^| shows the polymer length and ring size dis- 
tributions. In this simulation, 84.2% of the silane frag- 
ments form linear structures; the remainder form rings. 
For this system, the Si-R and the Si-Si bond lengths are 
2.09 ± 0.07 A and 2.32 ± 0.08 A. The R-Si-R and Si-Si- 
Si bond angles are 107.2 ± 5.0 and 111.1 ± 6.7 degrees, 
respectively. 




FIG. 2. Resulting polysilane structure at density 2.0 ■ 10 4 
monomers per A 3 . For clarity, the alkyl groups are not shown. 



FIG. 4. Resulting polysilane structure at density 5.0 ■ 10 
monomers per A 3 . For clarity, the alkyl groups are not shown. 
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FIG. 3. Histograms showing polymer length (a) and 
ring size distribution (b) for polysilane consisting of 1000 
monomers at density 2.0 ■ 10~ 4 fragments per A 3 . 

To study the effect of the monomer density, an addi- 
tional polysilane is simulated at an increased density of 
5.0 • 10 -3 monomers per A 3 . This network is shown in 
Fig. ||; 92.3% of all Si-Si bonds are formed. In this case, 






- 7 




12 




(b) " 


- 6 




10 








- 5 


'ence 


8 








- 4 


:cun 


6 








- 3 


o 


4 








- 2 

- 1 




2 






I on" 


- 







, n , 







5 10 15 20 25 30 35 40 
polymer length 



2 4 6 8 10 12 
ring size 



FIG. 5. Histograms showing polymer length (a) and ring 
size statistics (b) for polysilane consisting of 1000 monomers 
at density 5.0 • 10 -3 monomers per A 3 . 



B. Polysilyne 

We also performed one simulation involving 1000 SiR 
monomers, no SiR2 monomers being present at a density 



4 



of 2.0 • 10~ 4 monomer fragments per A 3 . The final net- 
work is shown in Fig. ^; 77.9% of all possible Si-Si bonds 
were formed. Since silyne is trifunctional, the genus of 
the structures can now reach large values. The final net- 
work consists of 96 separate clusters, ranging in size from 
two up to 29 monomers. Fig. shows a histogram of the 
cluster sizes. Of these clusters, only three are chains, 
31 possess a monocyclic structure and the remaining are 
structures with g > 1. Fig [?]b shows a histogram of the 
genus numbers. The structure is not dendritic: by far the 
largest fraction of monomers are part of rings; the frac- 
tion of monomers in chains or side-chains is 5.6% only. 
Moreover, the chains and side-chains are short, consist- 
ing at most of three monomers. For this network, the 
Si-R and the Si-Si bond lengths are 2.09 ± 0.08 A and 
2.33 ± 0.09 A, respectively. The Si-Si-Si bond angle is 
106.7 ± 10.0 degrees. 




FIG. 6. Polysilyne network at a density of 2.0 ■ 10 4 
monomers per A 3 . For clarity, the alkyl groups are not shown. 
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FIG. 7. Histogram of cluster size (a) and genus number (b) 
of a silyne network consisting of 1000 monomers at density 
2.0 ■ 10 -4 monomers per A 3 . 



To study the effect of the monomer density, an addi- 
tional silyne network is generated at an increased density 
of 8.0 • 10~ 3 monomers per A 3 . This network is shown 
in Fig. |[ In this case, 88.8% of all possible Si-Si bonds 
are formed. All monomers form one large cluster with a 
genus of 333. The structure consists mostly of connected 
rings; 94.4% of the monomers are part of rings. The ring 
statistics are shown in Fig. ^| 




FIG. 8. Polysilyne network at a density of 8.0 ■ 10 
monomers per A 3 . For clarity, the alkyl groups are not shown. 




FIG. 9. Ring size histogram for a polysilyne cluster con- 
sisting of 1000 monomers. This cluster was generated at a 
density of 8.0 ■ 10 -3 monomers per A 3 . 



C. Hybrid silane/silyne networks 

We also performed simulations of hybrid silane/silyne 
networks, either rich in silylene units (containing 800 
SiR2 and 200 SiR monomers) or rich in branching points 
(containing 200 S1R2 and 800 SiR monomers). The den- 
sity was 2.0 • 10 -4 monomers per A 3 . For both net- 
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works, approximately 80% of the possible Si-Si bonds are 
formed. The SiR,2-rich network consists of 175 separate 
clusters, ranging in size from one up to 17 monomers. 
In Fig. |ic|a, a histogram of the cluster sizes is shown. 
Of these clusters, 60 are linear structures, 98 are rings, 
the rest are structures with g > 1. See Fig. [To| b for 
a histogram of the genus numbers found. For this net- 
work, the typical chain length ranges from two to four 
monomers; rings consist typically of five to six silicon 
atoms. 

The silyne-rich network consists of 108 separate clus- 
ters, ranging in size from one up to 37 monomers. In 
Fig. [Tl|a, we show a histogram of the cluster sizes found. 
Of these clusters, 17 are linear structures, 31 are rings, 
the rest are structures with g > 1. See Fig. [Tl| b for a his- 
togram of the genus numbers found. For this network, 
the typical chain length is two monomers; rings typically 
contain five to six silicon atoms. 

Both networks show structures that are not dendritic. 
For the silane-rich network, 31.8% of the monomers form 
chains. For the silyne-rich network, this fraction is 11.6%. 
Most of the monomers are incorporated in rings. 
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FIG. 10. Histogram of cluster size (a) and genus number 
(b) for a silane-rich network, consisting of 800 silane and 200 
silyne monomers. 
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FIG. 11. Histogram of cluster size (a) and genus number 
(b) for a silyne-rich network, consisting of 200 silane and 800 
silyne monomers. 



V. DISCUSSION AND CONCLUSION 

The Wurtz-type polymerization of dichlorodialkylsi- 
lanes and trichloroalkylsilanes is a complex process . 
The polymerization proceeds by a chain growth process 
at the alkali metal surface, which plays an essential role 
in the formation of polymers of high molecular weight. 
This is for instance indicated by the fact that employ- 
ment of homogeneous coupling agents leads to formation 
of only oligomeric materials p|,p0||. Which products ac- 
tually are obtained also depends on parameters as the 
nature of the alkali metal, the solvent and the type of 
organic side chain. Moreover, after the initial stages of 
the formation of linear polymers, secondary reactions in- 
cluding depolymerization, backbiting by silyl anions and 
redistribution of chain lengths are known to occur . 
For polysilynes, secondary reactions such as ring open- 
ing reactions (by analogy with spirosilanes, see Ref. fbl), 
ring-size redistribution and coupling to linear moieties 
may be of interest. It is a virtually impossible task to 
perform a simulation of a Wurtz-type polymerization at 
a metal surface, hereby also taking all these contribut- 
ing factors into account. Nevertheless, as the approach 
we used led to results which adequately reflect experi- 
mental findings, we believe that our results are sound 
and provide insight into the structure and properties of 
polysilanes, polysilynes and hybrid polysilane/polysilyne 
networks. 

For polysilane, for instance, the simulations adequately 
predict the effect of the monomer density on the nature of 
the silicon compounds formed. At low density, the system 
has a strong tendency to form five- and six-membered 
rings, see Fig. [| This finding is in accordance with 
the experimental observation that formation of cyclic 
silicon-based compounds is favored by slow addition of 
the monomer to the alkali metal, i.e. by maintaining a low 
monomer concentration |21 2^]. The ring size distribu- 
tion shown in Fig. |3|, which indicates that five-membered 
rings are formed in excess, is in agreement with the ex- 
perimental distribution such as obtained upon reaction 
of a dialkyldichlorosilane under kinetic control It 
is noteworthy that in an equilibrium distribution, which 
arises from redistribution of the kinetic mixture upon 
use of excess alkali metal, the six-membcrcd ring dom- 
inates. As expected, the amount of long linear frag- 
ments increases strongly when the monomer density is 
increased, see Fig. 0. At this higher concentration five- 
and six-membered rings are formed in equal amounts. 
Furthermore, the length of the linear fragments is rather 
moderate, which may support the idea that the presence 
of an alkali metal surface is essential for the formation of 
high molecular weight polymers [pppj] . 

For the polysilyne, the simulations also predict the 
formation of cyclic structures, see Figs. |^-^[ At low 
concentration virtually all monomers are incorporated 
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into rings. Notwithstanding that monocyclic systems are 
formed in the largest amount, the number of mutually 
bonded rings is already considerable. At high monomer 
density, a very complex network of genus 333 is obtained. 
Of the monomers, 94.4% is incorporated in rings, which 
mainly contain five or six silicon atoms. Thus, whereas 
in the linear case a higher concentration leads to a larger 
amount of linear structures, for the polysilynes the sili- 
con atoms are still virtually exclusively found in cyclic 
structures. A few linear fragments are however also 
present. The simulated structure is in close agreement 
with the characteristics of the material that was isolated 
after the first stage of a polymerization of trichloroalkyl- 
silanes [^3|. Our simulations indicate that the structure 
of polysilynes can be viewed as a network of silicon atoms 
comprising fused cyclic structures and cyclic structures 
mutually connected by a single Si-Si bond. There is 
neither indication of formation of a regular sheetlike ar- 
rangement of silicon atoms nor of a dendritic struc- 
ture |p"3f , at least in our simulations where there is no 
difference in reactivity between Si-SiCl-Si groups and ter- 
minal SiCl groups. 

The results obtained for the hybrid polysilane / polysi- 
lyne networks do not substantially differ from those ob- 
tained for the linear or fully branched systems. With 20% 
SiR monomers present, mainly short linear and mono- 
cyclic compounds are formed, see Fig. [It]. However, the 
average cluster size is larger than in the linear case, and 
already some clusters of genus two and three, i.e. bi- 
cyclic and tricyclic systems are found. When the num- 
ber of branching points is increased to 80%, the average 
cluster size and the genus increase considerably. This 
demonstrates once again that incorporation of trifunc- 
tional (SiR) monomers leads to a structure with inter- 
connected and fused rings. Hence, the results for the 
hybrid polysilane and polysilyne networks show that this 
class of materials is built up from cyclic structures inter- 
connected via short linear chains. 

An interesting question is to what extent the structure 
of the polysilyne depicted in Fig. |^ is consistent with the 
flexibility of polysilynes such as indicated by its thermal 
properties, i.e. whether the structure allows for confor- 
mational changes affecting the degree of cr-conjugation. 
Although at first sight the network of fused cyclic struc- 
tures seems to be quite rigid, there are a number of fea- 
tures which suggest a certain degree of flexibility. Firstly, 
some (short) linear fragments and extended cyclic struc- 
tures, both of which imply conformational flexibility, are 
distinguishable in Fig. ^|. Secondly, there is also a num- 
ber of rings which are interconnected by a single Si-Si 
bond. This should give the possibility for rings to rotate 
with respect to each other; i.e., the structure depicted in 
Fig. ^| is not rigid. A third factor which may contribute 
to the flexibility of the networks is that silicon rings are 
known to be highly flexible. In this context it is of interest 
that even in the solid state dodecamethylcyclohexasilane 



SigMei2 undergoes a rapid ring inversion p4| . While it 
is not likely that ring inversions occur for fused ring sys- 
tems, they might be possible for rings connected to linear 
chains or connected to other rings by a single Si-Si bond. 

Hence, the networks obtained in our simulations are 
already flexible. However, there are a number of reasons 
why experimentally prepared polysilynes can even more 
easily undergo conformational changes. A preferred poly- 
mer extension at the termini of the already formed poly- 
mer, as discussed above, will enhance chain growth rather 
than branching. This will have the consequence that 
polysilynes in reality may contain more linear chains and 
larger cyclic structures than depicted in Fig. ^. Another 
factor is the size of the organic side group R. While in 
the simulations R has the dimensions of a methyl group, 
much more bulky side groups such as hexyl, isobutyl and 
phenyl substituents have been used in experimental stud- 
ies HJtJ ■ If bulky side groups are present the network is 
expected to be less dense and to incorporate more linear 
fragments than in the present case. 

A question left to be answered is whether polysilynes 
obtained by a Wurtz-type condensation should be re- 
garded as one, two or three-dimensional silicon materials. 
At first sight, the silicon network in Fig. |^ percolates in 
three dimensions. However, the description of the sim- 
ulated polysilyne and hybrid polysilane/polysilyne net- 
works by g = 1 + e — n, which is Euler's equation for 
two-dimensional networks ]25[ | , implies that topologically 
the silicon backbones of these materials are best regarded 
as two-dimensional systems. This is consistent with the 
electronic spectra which approach that of an indirect 
band gap semiconductor j6[j|]. 

In summary, we have shown that the formation 
and structure of polysilanes, polysilynes and hybrid- 
polysilane/polysilyne networks can be adequately de- 
scribed by a simulation based on rare event dynamics 
of diffusion and reaction steps. The simulations indicate 
that ring formation is an important factor for all three 
types of silicon materials. Insight has been obtained in 
the conformational flexibility of polysilynes and polysi- 
lane/polysilyne materials such as inferred from experi- 
mental studies. It is anticipated that simulations can pro- 
vide an even more comprehensive picture when variation 
of the size of the organic substituent R and differences 
in reactivity of different types of Si-Cl functionalities are 
implemented. This would be the subject of further work. 
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